A cost-effective blood DNA methylation-based age estimation method in domestic cats, Tsushima leopard cats (Prionailurus bengalensis euptilurus) and Panthera species, using targeted bisulphite sequencing and machine learning models.

Individual age can be used to design more efficient and suitable management plans in both in situ and ex situ conservation programmes for targeted wildlife species. DNA methylation is a promising marker of epigenetic ageing that can accurately estimate age from small amounts of biological material, which can be collected in a minimally invasive manner. In this study, we sequenced five targeted genetic regions and used 8-23 selected CpG sites to build age estimation models using machine learning methods at only about $3-7 per sample. Blood samples of seven Felidae species were used, ranging from small to big, and domestic to endangered species: domestic cats (Felis catus, 139 samples), Tsushima leopard cats (Prionailurus bengalensis euptilurus, 84 samples) and five Panthera species (96 samples). The models achieved satisfactory accuracy, with the mean absolute error of the most accurate models recorded at 1.966, 1.348 and 1.552 years in domestic cats, Tsushima leopard cats and Panthera spp. respectively. We developed the models in domestic cats and Tsushima leopard cats, which were applicable to individuals regardless of health conditions; therefore, these models are applicable to samples collected from individuals with diverse characteristics, which is often the case in conservation. We also showed the possibility of developing universal age estimation models for the five Panthera spp. using only two of the five genetic regions. We do not recommend building a common age estimation model for all the target species using our markers, because of the degraded performance of models that included all species.

Non-invasive age estimation based on faecal DNA using methylation-sensitive high-resolution melting for Indo-Pacific bottlenose dolphins.

Age is necessary information for the study of life history of wild animals. A general method to estimate the age of odontocetes is counting dental growth layer groups (GLGs). However, this method is highly invasive as it requires the capture and handling of individuals to collect their teeth. Recently, the development of DNA-based age estimation methods has been actively studied as an alternative to such invasive methods, of which many have relied on used biopsy samples. However, if DNA-based age estimation can be developed from faecal samples, age estimation can be performed entirely non-invasively. We developed an age estimation model using the methylation rate of two gene regions, GRIA2 and CDKN2A, measured through methylation-sensitive high-resolution melting (MS-HRM) from faecal samples of wild Indo-Pacific bottlenose dolphins (Tursiops aduncus). The age of individuals was known through conducting longitudinal individual identification surveys underwater. Methylation rates were quantified from 36 samples collected from 30 individuals. Both gene regions showed a significant correlation between age and methylation rate. The age estimation model was constructed based on the methylation rates of both genes which achieved sufficient accuracy (after LOOCV: MAE = 5.08, R2 = 0.33) for the ecological studies of the Indo-Pacific bottlenose dolphins, with a lifespan of 40-50 years. This is the first study to report the use of non-invasive faecal samples to estimate the age of marine mammals.

Characterization of Common Minke Whale (Balaenoptera Acutorostrata) Cell Lines Immortalized with the Expression of Cell Cycle Regulators.

Primary cultured cells cannot proliferate infinite. The overcoming of this limit can be classified as immortalization. Bypass of p16 senescence protein induces efficient immortalization various types of mammalians is previously reported. However, the Cetacea species is not known. Here, that common minke whale-derived cells can be immortalized with a combination of human genes, mutant cyclin-dependent kinase 4 (CDK4R24C ), cyclin D1, and Telomerase Reverse Transcriptase (TERT) is reported. These results indicate that the function of cell cycle regulators in premature senescence is evolutionarily conserved. This study describes the conserved roles of cell cycle regulators in the immortalization of cells from humans to Cetacea species. Furthermore, using RNA-seq based on next-generation sequencing, the gene expression profiles of immortalized cells are compared with parental cells as well as those immortalized with SV40 large T antigen, which is once a popular method for cellular immortalization. The profiling results show that newly established common minke-whale-derived immortaliozed cells have completely different profiles from SV40 cells. This result indicates that the expression of mutant CDK4, cyclin D1, and TERT enables to establish immortalized cell lines with different biological nature from SV40 expressing cells.

Age estimation of captive Asian elephants (Elephas maximus) based on DNA methylation: An exploratory analysis using methylation-sensitive high-resolution melting (MS-HRM).

Age is an important parameter for bettering the understanding of biodemographic trends-development, survival, reproduction and environmental effects-critical for conservation. However, current age estimation methods are challenging to apply to many species, and no standardised technique has been adopted yet. This study examined the potential use of methylation-sensitive high-resolution melting (MS-HRM), a labour-, time-, and cost-effective method to estimate chronological age from DNA methylation in Asian elephants (Elephas maximus). The objective of this study was to investigate the accuracy and validation of MS-HRM use for age determination in long-lived species, such as Asian elephants. The average lifespan of Asian elephants is between 50-70 years but some have been known to survive for more than 80 years. DNA was extracted from 53 blood samples of captive Asian elephants across 11 zoos in Japan, with known ages ranging from a few months to 65 years. Methylation rates of two candidate age-related epigenetic genes, RALYL and TET2, were significantly correlated with chronological age. Finally, we established a linear, unisex age estimation model with a mean absolute error (MAE) of 7.36 years. This exploratory study suggests an avenue to further explore MS-HRM as an alternative method to estimate the chronological age of Asian elephants.

Association between arginine vasopressin receptor 1A (AVPR1A) polymorphism and inequity aversion.

Although numerous studies have focused on brain functions related to inequity aversion, few have examined its genetic basis. Here, we show the association between estimated inequity aversion and polymorphisms in three genes associated with human sociality. Non-student adult participants took part in five economic game experiments on different days. Disadvantageous inequity aversion (DIA) and advantageous inequity aversion (AIA) were calculated from behavioural responses using Bayesian estimation. We investigated the association between genetic polymorphisms in the oxytocin receptor (OXTR rs53576), arginine vasopressin receptor 1A (AVPR1A RS3) and opioid receptor mu 1 (OPRM1 rs1799971) and inequity aversion. Regarding AVPR1A RS3, participants with the SS genotype had higher AIA than those with the SL or LL genotypes, but no association was found for DIA. Moreover, we observed no aversion associations for OXTR rs53576 or OPRM1 rs1799971. The results suggest that AVPR1A plays an important role in aversion when one's own gain is greater than that of others. Our findings may provide a solid theoretical basis for future studies on the relationship between genetic polymorphisms and inequity aversion.

Genetic dissection of behavioral traits related to successful training of drug detection dogs.

Drug detection dogs play integral roles in society. However, the interplay between their behaviors and genetic characteristics underlying their performance remains uninvestigated. Herein, more than 120,000 genetic variants were evaluated in 326 German Shepherd or Labrador Retriever dogs to profile the genetic traits associated with various behavioral traits related to the successful training of drug detection dogs. Behavioral breed differences were observed in 'friendliness to humans' and 'tolerance to dogs.' A genome-wide association study within both breeds identified 11 regions potentially associated with drug detection dog characteristics as well as 'interest in the target' and 'friendliness to humans,' which are related to drug detection abilities. Among them, 63 protein-coding genes, including Atat1 and Pfn2 known to be associated with anxiety-related or exploration behavior in mice, respectively, were located surrounding the identified candidate polymorphisms. This study highlights genetic characteristics associated with behavioral traits that are important for the successful training of drug detection dogs. Thus, these findings may facilitate improved breeding and training of these dogs.

Age estimation based on blood DNA methylation levels in brown bears.

Age is an essential trait for understanding the ecology and management of wildlife. A conventional method of estimating age in wild animals is counting annuli formed in the cementum of teeth. This method has been used in bears despite some disadvantages, such as high invasiveness and the requirement for experienced observers. In this study, we established a novel age estimation method based on DNA methylation levels using blood collected from 49 brown bears of known ages living in both captivity and the wild. We performed bisulfite pyrosequencing and obtained methylation levels at 39 cytosine-phosphate-guanine (CpG) sites adjacent to 12 genes. The methylation levels of CpGs adjacent to four genes showed a significant correlation with age. The best model was based on DNA methylation levels at just four CpG sites adjacent to a single gene, SLC12A5, and it had high accuracy with a mean absolute error of 1.3 years and median absolute error of 1.0 year after leave-one-out cross-validation. This model represents the first epigenetic method of age estimation in brown bears, which provides benefits over tooth-based methods, including high accuracy, less invasiveness, and a simple procedure. Our model has the potential for application to other bear species, which will greatly improve ecological research, conservation, and management.

Immortalization of primary cells derived from the endangered Ryukyu long-furred rat.

The Ryukyu long-furred rat is an endangered species confined to the southernmost three small islands of Japan (Amami-Oshima, Tokunoshima, and Okinawa). Its population is rapidly decreasing because of roadkill, deforestation, and feral animals. To date, its genomic and biological information are poorly understood. In this study, we successfully immortalized Ryukyu long-furred rat cells by expressing a combination of cell cycle regulators, mutant cyclin-dependent kinase 4 (CDK4R24C) and cyclin D1, together with telomerase reverse transcriptase or an oncogenic protein, the Simian Virus large T antigen. The cell cycle distribution, telomerase enzymatic activity, and karyotype of these two immortalized cell lines were analyzed. The karyotype of the former cell line immortalized with cell cycle regulators and telomerase reverse transcriptase retained the nature of the primary cells, while that of the latter cell line immortalized with the Simian Virus large T antigen had many aberrant chromosomes. These immortalized cells would be valuable for studying the genomics and biology of Ryukyu long-furred rats.

Induced pluripotent stem cells of endangered avian species.

The number of endangered avian-related species increase in Japan recently. The application of new technologies, such as induced pluripotent stem cells (iPSCs), may contribute to the recovery of the decreasing numbers of endangered animals and conservation of genetic resources. We established novel iPSCs from three endangered avian species (Okinawa rail, Japanese ptarmigan, and Blakiston's fish owl) with seven reprogramming factors (M3O, Sox2, Klf4, c-Myc, Nanog, Lin28, and Klf2). The iPSCs are pluripotency markers and express pluripotency-related genes and differentiated into three germ layers in vivo and in vitro. These three endangered avian iPSCs displayed different cellular characteristics even though the same reprogramming factors use. Japanese ptarmigan-derived iPSCs have different biological characteristics from those observed in other avian-derived iPSCs. Japanese ptarmigan iPSCs contributed to chimeras part in chicken embryos. To the best of our knowledge, our findings provide the first evidence of the potential value of iPSCs as a resource for endangered avian species conservation.

Detection of anti-ebolavirus antibodies in Ghanaian pigs.

Some filoviruses such as ebolaviruses and marburgviruses, cause hemorrhagic fever in humans and nonhuman primates. Pigs are suggested to play a potential role in the filovirus ecology. We investigated the seroprevalence of filovirus infection in pigs in Ghana. Using a viral glycoprotein (GP)-based enzyme-linked immunosorbent assay, we detected filovirus-specific immunoglobulin G antibodies in 5 of 139 samples. These positive sera showed specificities to four different filovirus species. Particularly, two of the positive sera reacted to GPs of two African ebolaviruses (i.e., Ebola virus and Taï Forest virus) in Western blotting. Our results suggest that these Ghanaian pigs were exposed to multiple filoviruses and emphasize the importance of continuous monitoring of filovirus infection in pig populations in West African countries.

Cryopreservation Competence of Chicken Oocytes as a Model of Endangered Wild Birds: Effects of Storage Time and Temperature on the Ovarian Follicle Survival.

For the conservation of endangered avian species, developing gamete preservation technologies is essential. However, studies in oocytes have not been widely conducted. In this study, assuming that the ovaries are transported to a research facility after death, we investigated the effect of ovary storage on oocytes for the purpose of cryopreserving avian female gametes by using a chicken as a model of endangered avian species. After excision, the ovaries were stored at either a low temperature (4 °C) or room temperature for 1-3 days. Ovarian follicles stored under different conditions for each period were examined by neutral red staining, histology, and gene and protein expression analysis. In addition, the pH of the storage medium after preserving the ovaries was measured. Then, ovarian tissues were vitrified to determine the cryopreservation competence. Storing the ovarian tissues at 4 °C kept the follicles viable and morphologically normal for 3 days with slow decline. In contrast, although different storage temperature did not influence follicle viability and morphology after only 1 day of storage, ovarian tissues stored at room temperature rapidly declined in structurally normal follicles, and viable follicles were rarely seen after 3 days of storage. Gene and protein expression analysis showed that apoptosis had already started on the first day, as shown by the higher expression of CASP9 under room temperature conditions. Furthermore, high expression of SOD1 and a rapid decline of pH in the storage medium under room temperature storage suggested the influence of oxidative stress associated with low pH in this condition on the follicle survivability in hen ovarian tissues. Our cryopreservation study also showed that ovarian tissues stored at 4 °C could recover after cryopreservation even after 3 days of storage. The described storage conditions and cryopreservation methods, which preserve chicken follicle survival, will lay the foundation of ovarian tissue preservation to preserve the fertility of wild female birds.

Comparison of the Genetic Diversity of the Captive and Wild Populations of the Tsushima Leopard Cat Using a GRAS-Di Analysis.

The Tsushima leopard cat (Prionailurus bengalensis euptilurus) (TLC) is a regional population of the Amur leopard cat (P. bengalensis euptilurus) that lives only on the Tsushima Island in Japan and is threatened with extinction. Because the TLC population is small, genetic management is important. In this study, we obtained the draft genome of the TLC and identified single-nucleotide polymorphism (SNP) markers using a genotyping by random amplicon sequencing-direct (GRAS-Di) analysis. We genotyped 31 captive individuals and 50 wild individuals, of which 48 were from a previous study. The identified SNPs were used to clarify the genetic diversity and genetic structure of the wild and captive populations of the TLC. The size of the genome was estimated to be about 2.42 Gb. The number of SNP markers developed was 139, and although PID and probability of exclusion obtained using these SNP markers were not as high as those reported in the studies of other wild species, these SNP markers could be used to identify individuals and parentage. Moreover, the genetic diversity indices of the captive population were similar to those of the wild population. These SNP markers will be useful for understanding the ecology of the TLC and planning conservation strategies.

Chlamydial species among wild birds and livestock in the foothills of Mt. Afadjato, Ghana.

The members of family Chlamydiaceae have a broad host range and cause many kinds of diseases in humans and animals. Several cases of Chlamydiaceae being detected in atypical hosts have been reported recently. Consequently, cross-species monitoring of Chlamydia in wildlife and livestock is pertinent for public health, animal hygiene and wildlife conservation. In this study, we conducted molecular surveillance of Chlamydia in wild birds and livestock around a small village in the foothills of Mt. Afadjato, Ghana where direct contact between wildlife and livestock occurs. Among 29 captured wild birds and 63 livestock, 5 sheep, 30 goats and 28 chickens, the positive ratios of Chlamydia were 24.1%, 40.0%, 43.3% and 26.9%, respectively. Chlamydia pecorum was detected in wild birds, goats, sheep and chickens. On the basis of the variable domain 2 region of ompA, several samples from different hosts showed identical sequences and were phylogenetically located to the same clusters. In addition, using ompA, C. psittaci, C. abortus and C. gallinacea were also detected in this small habitat. Further genetic and pathogenic analyses of the chlamydial distribution in this area, which represents the interface of wild and domestic animal interactions, may improve our knowledge of their transmission among different hosts.

Association Between Polymorphism (5-HTTLPR) of the Serotonin Transporter Gene and Behavioral Response to Unfair Distribution.

Behavioral responses to unfair distribution have been measured mainly using the Ultimatum Game (UG). Recent studies examining the biological basis of behavioral responses to unfair distribution have focused attention on the role of the serotonin transporter gene. However, studies, to date, have been conducted on non-Asians, and it has not been confirmed whether similar results can be seen in other ethnic groups. It has also been shown that behavioral responses to unfair distribution are not only seen in the case of victims themselves but also in the case of third parties not directly affected. This study aimed to determine whether the results of the previous study would be replicated in an Asian population and whether the serotonin transporter gene would also be associated with behavior toward unfair distribution by third parties. We examined the association between polymorphisms (5-HTTLPR) of the serotonin transporter gene and participants' behavior in the UG and the third-party punishment game (TPPG). The results did not show an association between punishment for unfair proposals in the TPPG and genetic polymorphisms, while participants with the SL/LL genotype were more likely to reject unfair offers in the UG than those with the SS genotype. These results indicate that 5-HTTLPR is associated with behavior when unfair intentions are directed at oneself.

Generation and Gene Expression Profiles of Grevy's Zebra Induced Pluripotent Stem Cells.

Induced pluripotent stem cells (iPSCs) can serve as a biological resource for functional and conservation research for various species. This realization has led to the generation of iPSCs from many species, including those identified as endangered. However, the understanding of species variation in mammalian iPSCs remains largely unknown. To gain insight into species variation in iPSCs, we generated iPSCs from a new species Grevy's zebra (Equus grevyi; gz-iPSCs), which has been listed as endangered in the IUCN (International Union for Conservation of Nature) Red List. We isolated primary fibroblast cells from an individual and successfully reprogrammed them into iPSCs. The generated gz-iPSCs continued to grow under primed-type culture condition and showed pluripotency and differentiation potential. To describe the molecular characteristics of gz-iPSCs, we performed RNA sequencing analysis. The gz-iPSC transcriptome showed robust expression of pluripotency-associated genes reported in human and mouse, suggesting evolutionary conservation among the species. This study provides insight into the iPSCs from a rare species and helps the understanding of the gene expression basis underlying mammalian pluripotent stem cells.

Genetic and Phenotypic Characterization of Domestic Geese (Anser anser) in Egypt.

The objectives of this study were to achieve phenotypic characterization of three domestic Egyptian goose populations collected from three different geographical zones (Kafr El-Sheikh, Fayoum and Luxor) and to perform genetic characterization of these three populations based on mtDNA D-loop and 12 microsatellite markers. The body measurements of 402 domestic mature geese belonging to these three governorates showed that the lengths of the head, culmen and tarsus and the live body weight varied significantly among the three studied Egyptian goose populations. After alignment of a 710-base-pair segment of the goose mtDNA control region, there was a single haplotype in the three Egyptian goose populations, indicating the same maternal origins. The genotyping of the 12 microsatellite markers showed low diversity indices, including average observed (NA) and effective (NE) number of alleles and observed (HO) and expected heterozygosity (HE) (3.333, 1.760, 0.277 and 0.352, respectively), and a high inbreeding coefficient (FIS = 0.203) across the three Egyptian goose populations. The high inbreeding and low genetic and morphological differentiation of Egyptian geese could be corrected by establishing a large base population through capturing small populations with the highest genetic variation. The findings of the current study can therefore serve as an initial guide to design further investigations for developing conservation programs of Egyptian geese genetic resources.

Age estimation using methylation-sensitive high-resolution melting (MS-HRM) in both healthy felines and those with chronic kidney disease.

Age is an important ecological tool in wildlife conservation. However, it is difficult to estimate in most animals, including felines-most of whom are endangered. Here, we developed the first DNA methylation-based age-estimation technique-as an alternative to current age-estimation methods-for two feline species that share a relatively long genetic distance with each other: domestic cat (Felis catus; 79 blood samples) and an endangered Panthera, the snow leopard (Panthera uncia; 11 blood samples). We measured the methylation rates of two gene regions using methylation-sensitive high-resolution melting (MS-HRM). Domestic cat age was estimated with a mean absolute deviation (MAD) of 3.83 years. Health conditions influenced accuracy of the model. Specifically, the models built on cats with chronic kidney disease (CKD) had lower accuracy than those built on healthy cats. The snow leopard-specific model (i.e. the model that resets the model settings for snow leopards) had a better accuracy (MAD = 2.10 years) than that obtained on using the domestic cat model directly. This implies that our markers could be utilised across species, although changing the model settings when targeting different species could lead to better estimation accuracy. The snow leopard-specific model also successfully distinguished between sexually immature and mature individuals.

Personality, subjective well-being, and the serotonin 1a receptor gene in common marmosets (Callithrix jacchus).

Studies of personality traits in common marmosets (Callithrix jacchus) indicate that there are five or six constructs-Sociability, Dominance, Neuroticism, Openness, and two related to Conscientiousness. The present study attempted to determine whether our earlier study of laboratory-housed individuals only yielded three-Dominance, Sociability, and Neuroticism-because of a low amount of between-subjects variance. To do so, we increased our sample size from 77 to 128. In addition, we ascertained the reliability and validity of ratings and whether polymorphisms related to the serotonin 1a receptor were associated with personality. We found Sociability, Dominance, and Negative Affect factors that resembled three domains found in previous studies, including ours. We also found an Openness and Impulsiveness factor, the latter of which bore some resemblance to Conscientiousness, and two higher-order factors, Pro-sociality and Boldness. In further analyses, we could not exclude the possibility that Pro-sociality and Boldness represented a higher-level of personality organization. Correlations between personality factors and well-being were consistent with the definitions of the factors. There were no significant associations between personality and genotype. These results suggest that common marmoset personality structure varies as a function of rearing or housing variables that have not yet been investigated systematically.

Genetic diversity and population structure in the Ryukyu flying fox inferred from remote sampling in the Yaeyama archipelago.

The Ryukyu flying fox (Pteropus dasymallus) is distributed throughout the island chain spanning across southern Japan, Taiwan, and possibly the Philippines. Although P. dasymallus is listed as VU (vulnerable) in the IUCN Red List, only few genetic works have been conducted to support its conservation. In this study we analyzed 19 markers (mtDNA haplotypes and 18 microsatellite markers) to evaluate genetic diversity and investigate the genetic structure of this species. mtDNA analysis was conducted with 142 DNA remote samples, mostly from faeces, and wing tissues collected on eight islands (Miyako, Ishigaki, Kohama, Kuroshima, Hateruma, Taketomi, Iriomote, Yonaguni). 39 haplotypes were identified in 526bp of the control region, and haplotype network showed no clear genetic structure. Microsatellite analysis was also conducted with 155 samples collected on six islands (Miyako, Ishigaki, Kohama, Taketomi, Iriomote, Yonaguni). It showed that the Yonaguni population exhibits low genetic diversity, high inbreeding, and clear genetic differentiation from other populations. Gene flow between Ishigaki and Miyako through small stepstone islands might be preventing inbreeding of the Miyako population. We provide for the first time indirect proof of long-distance inter-island dispersal in the Ryukyu flying fox and revealed genetic diversity, gene flow and genetic differentiation among the archipelago's populations. These results will be useful for delineating conservation units and designing specific conservation policies for each island based on metapopulation genetic structure.